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Abstract 

We study large-scale dynamo action due to turbulence in the presence of a linear shear flow. Our treat- 
ment is quasilinear and equivalent to the standard 'first order smoothing approximation'. However it is non 
perturbative in the shear strength. We first derive an integro-differential equation for the evolution of the 
mean magnetic field, by systematic use of the shearing coordinate transformation and the Galilean invari- 
ance of the linear shear flow. We show that, for non helical turbulence, the time evolution of the cross-shear 
components of the mean field do not depend on any other components excepting themselves; this is valid 
for any Galilean-invariant velocity field, independent of its dynamics. Hence, to all orders in the shear pa- 
rameter, there is no shear-current type effect for non helical turbulence in a linear shear flow, in quasilinear 
theory in the limit of zero resistivity. We then develop a systematic approximation of the integro-differential 
equation for the case when the mean magnetic field varies slowly compared to the turbulence correlation 
time. For non-helical turbulence, the resulting partial differential equations can again be solved by making 
a shearing coordinate transformation in Fourier space. The resulting solutions are in the form of shearing 
waves, labeled by the wavenumber in the sheared coordinates. These shearing waves can grow at early and 
intermediate times but are expected to decay in the long time limit. 
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I. INTRODUCTION 



The origin of large-scale magnetic fields in astrophysical systems from stars to galaxies is 
an issue of considerable interest. The standard paradigm involves dynamo amplification of seed 
magnetic fields due to turbulent flows which have helicity combined with shear. Shear flows and 
turbulence are ubiquitous in astrophysical systems although the turbulence in general may not 
be helical. However the presence of shear by itself may open new pathways to the operation of 
large-scale dynamos, even if the turbulence lacks a coherent helicity yJ-lSD. The evidence for such 
large-scale dynamo action under the combined action of non helical turbulence and background 
shear flow comes mainly from several direct numerical simulations yl |21]. How such a dynamo 
works is not yet clear. One possibility is the shear-current effect in which extra components 
of the mean electromotive force (EMF) arise due to shear, which couple components of the mean 
magnetic field parallel and perpendicular to the shear flow. However there is no convergence yet 
on whether the sign of the relevant coupling term is such as to obtain a dynamo; some analytic 
calculations ^,2] and numerical experiments [I] find that the sign of the shear-current term is 
unfavorable for dynamo action. 

n 

In an earlier paper [8J (Paper I), we had outlined briefly a quasilinear theory of dynamo action 
in a linear shear flow of an incompressible fluid which has random velocity fluctuations due either 
to freely decaying turbulence or generated through external forcing. Our analysis did not put any 
restrictions on the strength of the shear, unlike earlier analytic work which treated shear as a small 
perturbation. We arrived at an integro-differential equation for the evolution of the mean magnetic 
field and argued that the shear-current assisted dynamo is essentially absent in quasilinear theory 
in the limit of zero resistivity. In the present paper we give detailed derivations of the main results 
of Paper I. We also extend our work further by deriving differential equations for the mean field, 
in the limit when the correlation time of the turbulence is much smaller than the time-scale over 
which the mean field varies. This allows us to solve for the mean field evolution in terms of 
the velocity correlation functions. We can draw some general conclusions on the shear dynamo 
independent of the exact velocity dynamics. In particular we note that the shear dynamo can lead 
to transient growth of large-scale fields in the form of shearing waves, but these waves ultimately 
decay, even in the absence of microscopic diffusion. 

In the section II we formulate the shear dynamo problem. Our theory is 'local' in character: In 
the laboratory frame we consider a background shear flow whose velocity is unidirectional (along 
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the X2 axis) and varies linearly in an orthogonal direction (the Xi axis). Section III outlines a 
quasilinear theory of the shear dynamo. Systematic use of the shearing transformation allows us 
to develop a theory that is non perturbative in the strength of the background shear. However, 
we ignore the complications associated with nonlinear interactions, hence magnetohydrodynamic 
(MHD) turbulence and the small-scale dynamo; so our theory is quasilinear in nature, equivalent 
to the 'first order smoothing approximation' (FOSA) Jq, 10]. The linear shear flow has a basic 



symmetry relating to measurements made by a special subset of all observers, who may be called 
comoving observers. This symmetry is the invariance of the equations with respect to a group 
of transformations that is a subgroup of the full Galilean group. It may be referred to as 'shear- 
restricted Galilean invariance', or Galilean invariance (GI). It should be noted that the laboratory 
frame and its set of comoving observers need not be inertial frames; in fact one of the main 
applications of GI is to the shearing sheet which is a rotating frame. We introduce and explore 
the consequences of GI velocity fluctuations in section IV. Such velocity fluctuations are not only 
compatible with the underlying symmetry of the problem, but they are expected to arise naturally. 
This has profound consequences for dynamo action, because the transport coefficients that define 
the mean EMF become spatially homogeneous in spite of the shear flow. The derivation of an 
integro-differential equation for the mean magnetic field is given in section V. We discuss a number 
of ways of approximating this equation in section VI, for slowly varying mean fields, all of which 
lead to the same set of partial differential equations for the mean-field. The mean field dynamics 
is further studied in section VII, and section VIII presents a discussion of the main results and the 
conclusions. 



n. THE SHEAR DYNAMO PROBLEM 

Let (ei, 62, 63) be the unit vectors of a Cartesian coordinate system in the laboratory frame, 
X = (Xi, X2, X3) the position vector, and r the time. The fluid velocity is given by (— 2^X162 + 
v), where A is the shear parameter (Oort's first constant) and v(X, r) is a randomly fluctuating 
velocity field. The total magnetic field, B'{X, r), obeys the induction equation: 

- 2AX,^^ B' + 2AB[e2 = Vx {vxB') + r]V'B' (1) 
The shear dynamo problem may be stated as follows: given some statistics of velocity fluctuations. 
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what can be said about the magnetic field? More specific questions may be posed: does the 
combined action of the background shear and random velocities lead to the growth of a large- 
scale component of the magnetic field (i.e. a turbulent dynamo)! In particular, is there turbulent 
dynamo action when the velocity fluctuations possess mirror-symmetry (i.e. when the velocity 
fluctuations are non helical)! 

A common approach to the problem is through the theory of mean-field electrodynamics . Here, 
the action of zero-mean velocity fluctuations {{v) = 0) on some seed magnetic field is assumed 
to produce a total magnetic field with a well-defined mean-field (B) and afluctuating-field (b): 

B' = B + b, {B') = B, (6) = (2) 

where ( ) denotes ensemble averaging in the sense of Reynolds. Applying Reynolds averaging to 
the induction equation ([T]), we obtain the following equations governing the dynamics of the mean 
and fluctuating magnetic fields: 



^ - 2AXi^ ] B + 2ABie2 = Vx£ + 7]V'B (3) 



2AXi-—]b + 2Abie2 = Vx(vxB) + Vx (vxb- £) + nV^b 

OT 0X2 J 

(4) 

where £ = (vxb) is, the mean EMF. The first step toward solving the problem is to calculate £ 
and obtain a closed equation for the mean-field, B(X, r). In the general case, it is necessary to 
specify the dynamics of v which could be influenced by Lorentz forces due to both B and b. 



III. QUASILINEAR THEORY 

To calculate the mean EMF we make some simplifying assumptions. We first make the quasi- 
linear approximation in solving equation dH) for b by dropping terms that that are quadratic in the 
fluctuations. Note that the dynamics of v is not prescribed; it does not imply absence of velocity 
dynamics. For instance, the fluid can be acted upon by Lorentz forces due to the magnetic field, 
Coriolis force as in the case of the shearing sheet, or buoyancy in a convective flow. In this paper 
we will not specify any particular dynamics for the velocity field. We also drop the resistive term 
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in the interests of simplicity of presentation. Setting 77 = may seem like a drastic step, but we 
would like to assure the reader that the theory can be reworked without this limitation and that our 
main conclusions carry through, even for r/ 7^ 0. In particular we recover the results of this paper 
in the limit r] 0. We note that the limit 77 — )• is also compatible with the physical situation 
in which the correlation times are small compared to the eddy turn-over timescale; so our theory 
is applicable when the 'first-order-smoothing-approximation' (FOSA) is valid. The fluctuating 
velocity field is assumed be incompressible (V- v = 0). This restriction is not crucial and may 
be lifted without much difficulty. 

The quasilinear approximation is equivalent to neglecting the effects of magnetohydrodynamic 
turbulence and small-scale dynamo action, for the determination of £. With these assumptions, 
the equation for b we will solve is 



- 2AX1-— 1 b + 2Ahe2 = Vx (vxB) 



[B-V)v - {vV)B (5) 



A. The shearing coordinate transformation 

Equation (|5]) is inhomogeneous in the coordinate Xi. It is convenient to exchange spatial inho- 
mogeneity for temporal inhomogeneity, so we get rid of the {Xid/dX2) term through a shearing 
transformation to new spacetime variables: 

Xi = Xi , X2=X2 + 2AtXi , X3 = X3 , t = T (6) 

Then partial derivatives transform as 



d -)- 2At ^ ^ + 2y4 (7) 

dXi dxi dx2 ' dX2 8x2 ' dX^ dxs ' dr dt 8x2 

We also define new variables, which are component-wise equal to the old variables: 

H{x,t) = B{X,t), h{x,t) = b{X,T), u{x,t) = v{X,t) (8) 

It is important to note that, just like the old variables, the new variables are expanded in the 
fixed Cartesian basis of the laboratory frame . For example, // = Hiei + H2e2 + H^e^, where 
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Hi{x,t) = Bi{X,T), and similarly for the other variables. In the new variables, equation ([5]) 
becomes, 

dh f d d \ f d d \ 

— + 2Ahie2 = {H- — + 2AtHi— ]u - [u- — + 2Atui— H (9) 

at \ ox 0x2 J \ ox 0x2 J 

Equation dH) for h{x, t) does not contain spatial derivatives of h, so it can be integrated directly. 
We are interested in the particular solution which vanishes at t = 0. The solutions for hi{x, t) and 
h-i{x, t) are: 



hi = [ dt'u[i[H'i+2At'6i2H[] - [ dt' [u'l + 2At'6i2u[] H[i (10) 
Jo Jo 



h3 = [ dt'u'^i[H[ + 2At'5i2H[] - [ dt' [u'i + 2At'6i2u[]H'^i (11) 
Jo Jo 

where primes denote evaluation at spacetime point {x,t'). We have also used notation Umi = 

(dum/dxi) and H,r,i = {dH„Jdxi). 

The equation for h2{x, t) involves hi{x, t); the solution is 



h2 = [ dt' u'^i[H[ + 2At'5i2H[] - [ dt' [u'l + 2At'6i2u[] H'^i - 2A [ dt' h[ (12) 
Jo Jo Jo 

We need to evaluate the integral 



t rt' 



Jo 



dt'h[ = I dt' I dt" H'i +2At"Si2H'l 



t rt' 

dt' / dt" 
Jo 



Ui + 2At Oi2U^ 



H 
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where the double-primes denote evaluation at spacetime point {x, t"). We now note that, for any 
function f{x, t), the double-time integral 



[ dt' [ dt" f{x,t") = [ dt" f{x,t") [ dt' 
Jo Jo Jo Jt" 

= fdt'{t-t')f{x,t') 
Jo 



dt" {t-t")f{x,t") 



reduces to a single-time integral, where in the last equality we have merely replaced the dummy 
integration variable t" by t'. Then 
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dt'h[ = / dt' {t-t')u[i[Hi + 2At'5i2H[] - / dt' {t - t') [u'l + 2At'6i2u[] H[i (14) 
Jo Jo 

can be used in equation (fT2l) to get an explicit solution for h2{x,t). Combining equations (flOl) . 

(fTTl) and (fT2l) we can write h{x,t) in component form as 

hm{x,t) = f dt' [u'^, - 2A{t-t')5^2u[i][H[ + 2At'5nH[] 
Jo 

- [ dt' [u[ + 2At'6i2u[][H'^i - 2A{t-t')6^2K] (15) 
Jo 

B. The mean EMF 

The expression in equation (fT5l) for h should be substituted in £ = (vxb) = (uxh). Follow- 
ing standard procedure, we allow ( ) to act only on the velocity variables but not the mean field; 
symbolically, it is assumed that (uuH) = (uu) H. Interchanging the dummy indices (/, m) in 
the last term of equation (fT5l) , the mean EMF is given in component form as 



[ dt' \Mx,t,t') - 2A{t-t')^u{x,t,t')][H'i + 2At'6i2H',] 
Jo 

/ dt' [fiimi{x,t,t') + 2At'6m2rjiii{x,t,t')][H'i^ - 2A{t-t')6i2H'^, 
Jo 



(16) 



where the transport coefficients, {a , 13 ,r]), are defined in terms of the uu velocity correlators by 



aii{x,t,t') = eijm{Uj{x,t)Uml{x,t')) 

^ii{x,t,t') = eij2{uj{x,t)uu{x,t')) 

r]imiix,t,t') = eiji{Uj{x,t)Um{x,t')) (17) 
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To obtain more specific expressions for the transport coefficients, we need to provide information 
on the uu velocity correlators. However, it is physically more transparent to consider velocity 
statistics in terms of vv velocity correlators, because this is referred to the laboratory frame instead 
of the sheared coordinates. By definition (eqn.[8]). 



Umix,t) = VmiX{x,t),t) (18) 



where 



Xi = xi, X2 = X2-2Atxi, X3 = X3, T = t (19) 
is the inverse of the shearing transformation given in equation Using 

d d d 

3^. = ax. - '^''''ox. 

the velocity gradient Umi can be written as 

f d d \ 

Uml = - 2Ar 5a —— Vm = Vml - 2At 611 Vm2 (21) 



ydXi 0X2/ 

where Vmi = (dvm/dXi). Then the transport coefficients are given in terms of the vv velocity 
correlators by 

%{x,tX) = t,,2[{v,{X,t)vii{X',t')) - 2At' 6n {v,{X,t)vi2{X',t'))] 
rfimi{x,t,t') = eiji{vj{X,t)vm{X',t')) (22) 

where X and X' are shorthand for 

X = {xi ,X2 — 2Atxi ,Xs) , X' = {xi ,X2 — 2At'xi ,X3) (23) 

Equation (fT6l) . together with (fTTI) or (|22l) . gives the mean EMF in general form. X and X' can 
be thought of as the coordinates of the origin, at times t and t' respectively, of an observer comoving 
with the background shear flow. Therefore the transport coefficients depend only on the velocity 
correlators measured by such an observer at the origin of her coordinate system. This fact will have 
profound consequences for dynamo action, when we consider G-invariant velocity correlators in 
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the next section. Before discussing the Galilean invariance of the linear shear flow, we derive the 
form of the mean EMF for a special case, when the velocity field is "delta-correlated-in-time". 



C. delta-correlated-in-time velocity correlator 

Although somewhat artificial, it is not uncommon to study dynamo action due to velocity fields 
whose correlation times are supposed so small that the two-point correlator taken between space- 
time points {R, r) and (R', r') is assumed to be 

{v,{R,t)v,{R',t')) = 6{r-T')T,,{R,R',T) (24) 
Incompressiblility implies that 

0; ^ = (25) 



We define 



The delta-function ensures that X and X' defined in equation (l23l) are equal to each other. Then 
the velocity correlators 



{v,{X,t)vj{X',t')) = 6{t-t')T,,{X,X,t) 

{v,{X,t)v,i{X',t')) = 6{t-t')T,,i{X,t) (27) 



Substitute equation (|27l) in equation (|22l) for the transport coefficients 



aii{x, t, t') = 5{t- t') eijm [Tjmi - 2At 5/i T^-^a] 
$ii{x,t,i:) = 6{t-t')e,j2[Tjii - 2At6nTju] 
rji^i{x,t,t') = 5{t-t')eijiT^m (28) 

and use these expressions in equation (fT6l ). The delta-function ensures that the integrals over time 
can all be performed explicitly, so the mean EMF is 
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8i = eijm[Tjmi - 2AtdnTj^2\[Hi + 2Atdi2H,] - eiji[Tjm + 2Atdm2Tji] Hi^ (29) 

It is useful to write the EMF in terms of the original variables and laboratory frame coordinates. 
To this end we transform 

Him = - 2ATdml^] Bi = Bim - 2At S^l B12 (30) 



where Bim — (dBi/dXm)- Then the explicit dependence of Ei on the shear parameter A cancels 
out, and the mean EMF assumes the simple form, 

^ijmTjml Bi €ijiTjm Bim (31) 

which is identical to the familiar expression in the absence of background shear. Therefore we 
conclude that, to obtain non trivial effects due to the shear flow, it is necessary to consider velocity 
correlators with non zero correlation times. Henceforth we shall consider the general case of finite 
velocity correlation times. 



IV. GALILEAN INVARIANCE 



The linear shear flow has a basic symmetry relating to measurements made by a special subset 
of all observers. We define a comoving observer as one whose velocity with respect to the labora- 
tory frame is equal to the velocity of the background shear flow, and whose Cartesian coordinate 
axes are aligned with those of the laboratory frame. A comoving observer can be labeled by the 
coordinates, $, = {^1,^2,^?,) with respect to the laboratory frame, of her origin at time r = 0. 
Different labels identify different comoving observers and vice versa. As the labels run over all 
possible values, they exhaust the set of all comoving observers. The origin of the coordinate axes 
of a comoving observer translates with uniform velocity; its position with respect to the origin of 
the laboratory frame is given by 

^c(t) = (6,6-2^x6,6) (32) 

An event with spacetime coordinates (X, r) in the laboratory frame has spacetime coordinates 
(JC , f) with respect to the comoving observer, given by 

X ^ X - Xc{t) , f ^ T -To (33) 
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where the arbitrary constant tq allows for translation in time as well. 



Let 



denote the total, the mean, the fluctuating mag- 



netic fields and the fluctuating velocity field, respectively, as measured by the comoving observer. 
These are all equal to the respective quantities measured in the laboratory frame: 



B\X,f),B{X,f)Mx,f),v{X,f)\ = [B\X,T),B{X,T),b{X,T),v{X,r)] (34) 

That this must be true may be understood as follows. Magnetic fields are invariant un- 
der non-relativistic boosts, so the total, mean and fluctuating magnetic fields must be the 
same in both frames. To see that the fluctuating velocity fields must be the same, we note 
that the total fluid velocity measured by the comoving observer is, by definition, equal to 
(^—2AXie2 + v{X, f) j . This must be equal to the difference between the velocity in the lab- 
oratory frame, {—2AXie2 + v(X,t)), and (— 2^4^162), which is the velocity of the comoving 
observer with respect to the laboratory frame. Using X = X —^1, we see that v{X, f) = v(X, r). 

The Galilean coordinate transformation given in equation (|33l) imphes that partial derivatives 
are related through 

9 d d d ^ d ^^^^ 

dX dX or Or 8X2 



Note that the combination (d/dr — 2AXid/dX2) = [d/df — 2AXid/dX2j is invariant in form. 
The other partial derivatives occurring in equations ©, ^ and (H)) are spatial derivatives which, 
by the second of equations (1351 ). are the same in both frames. Therefore equations ([T]), (|3]) and 
dH) are invariant under the simultaneous transformations given in equations (|33l) and (|34l ). We 
note that this symmetry property is actually invariance under a subset of the full ten-parameter 
Galilean group, parametrized by the five quantities (^1, ,^2, ^3, tq, for brevity we will refer to 
this restricted symmetry as Galilean invariance, or simply GI. 

There is a fundamental difference between the coordinate transformations associated with 
Galilean invariance (equation |33]) and the shearing transformation (equation |6l). The former re- 
lates different comoving observers, whereas the latter describes a time-dependent distortion of the 
coordinates axes of one observer. Comparing equation (|35]) with (|7]), we note that the relation- 
ship between old and new variables is homogeneous for the Galilean transformation, whereas it is 
inhomogeneous for the shearing transformation. 



11 



It is important to note that the laboratory frame and its set of comoving observers need not be 
inertial frames. Indeed, one of the main applications of our theory is to the shearing sheet which 
is a rotating frame providing a local description of a differentially rotating disc; in addition to 
other forces, the velocity field is affected by the Coriolis force. The only requirement is that the 
magnetic field satisfies the induction equation ([T]). 

A. Galilean-invariant velocity correlators 

Naturally occurring velocity fields are Galilean-invariant, and this has a strong impact on the 
velocity statistics. We consider the ra-point velocity correlator measured by the observer in the lab- 
oratory frame. Let this observer correlate Vj^ at spacetime location ri), with Vj^ at spacetime 
location T2), and so on upto Vj^ at spacetime location r„). Now consider a comoving 
observer, the position vector of whose origin is given by Xc{t) of equation (l32l) . An identical ex- 
periment performed by this observer must yield the same results, the measurements now made at 
the spacetime points denoted by {Ri + Xc{ti), ti) ; {R2 + X 

If the velocity statistics is GI, the n-point velocity correlator must satisfy the condition 



{Vj^{Rl,Ti) . . .Vj^{Rn,Tn)) = {Vj^{Ri + Xc{Ti),Ti)...Vj^{Rn + Xc{rn),Tn)) (36) 



for all {Ri, . . . Rn ; ri, . . . r„ ; ^). In quasilinear theory we require only the two-point velocity 
correlators, for which 



for all (i?, R! , r, r', ^). We also need to work out the correlation between velocities and their 
gradients: 



{v,{R,T)v,{R',r')) 



{viiR + X,(r), r) v,iR' + X,{r'), r')) 



(37) 



{viiR, T)vjiiR',r')) 



d 



{v,{R,t)v,{R',t')) 



dR' 



d 



{v,{R + X,(r), r) Vj{R! + X,(r'), r')) 



OR' 



{v,iR + Xe(r), r) v.iiR' + X,(r'), r')) 



(38) 
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If we now set 



R = R' = 0, T = t, t' = t' , (^1,6,6) = {xuX2,x^) (39) 
we will have 

Xe(r) = (a;i,a;2 -2Atxi,X3) , X^{t') = {xi,X2~2At'xi,x^) (40) 

Comparing equation (|40l) with equation (|23] ). we see that Xc(r) and Xc(r') are equal to X and 
X' , which are quantities that enter as arguments in the velocity correlators of equations (|22)) 
defining the transport coefficients. Hence, reading equations (|37l) and (|38l) from right to left, 
the velocity correlators. 



{v,{X,t)v,{X',t')) = {v,{0,t)v,{Q,t')) = R.j{t,t') 
{viiX,t)v,i{X',t')) = {v,{0,t)v,i{0,t')) = S^,i{t,t') (41) 

are independent of space, and are given by the functions, Rij{t, t') and Siji{t, t'). Symmetry and 
incompressiblity imply that Rij{t, t') = Rji{t', t) and Sijj{t, t') = . Note that the turbulence will, 
in general, be affected by the background shear and the velocity correlators will not be isotropic. 
In particular, Rij (t, t') will not be proportional to the unit tensor, 6ij. 



B. Galilean-invariant mean EMF 



The transport coefficients are completely determined by the form of the velocity correlator. 
Using equations (|4T]) in equations (|22b . we can see that the GI transport coefficients. 



0^il{t,t') — eijm[Sjmlit,t') — 2At' 5ii Sjm2{tTt')] 

Mt,t') = eij2[Sjii{t,t') - 2At'6nSju{t,t')] 
Vimi{t,t') = eijiR,m{t,t') (42) 

are independent of space. Galilean invariance is the fundamental reason that the velocity corre- 
lators, hence the transport coefficients, are independent of space. The derivation given above is 
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purely mathematical, relying on the basic freedom of choice of parameters made in equation (|39l ). 
but we can also understand the results more physically. X and X', as given by equation (l23l) . can 
be thought of as the location of the origin of a comoving observer at times t and t', respectively. 
Thus when the observer correlates velocities at X = Xc{t) and X' = Xc{t'), it will be the same 
as correlating the velocities at her origin, but at different times. Then GI implies that the velocity 
correlators must be equal to those measured by any comoving observer at her origin at times t 
and t' . In particular, this must be true for the observer in the laboratory frame, which explains 
equations (|4T]) . consequently equations (|42|) . 

We can derive an expression for the G-invariant mean EMF by using equations (|42l) for the 
transport coefficients in equation (fT6l) . The integrands can be simplified as follows: 

aii{t,t') [H[ + 2At'5i2H[] 
Pa{t,t') [H'i+2At'6i2H[] 

[ffimiit, t') + 2 At' 5^2 miiit, t')] = eiji [Rjmit, t') + 2At'5m2 Rjiit, t')] 

[%^2{t, t') + 2 At' 5^2 %i2{t, t')] H'^^ = ei^2 Sii [R^mit, t') + 2 At' 5,^2 Rji{t, t')] H'^^ 
Define 

Cjml{t,t') = Sjml{t,t') — 2A{t — t')Sm2 Sjil{t,t') 

D,^{t,t') = RjUt,t') + 2At'Sm2Rji{t,t') (43) 
The mean EMF can now be written compactly as 

S,{x,t) = ei,m [ dt' Cjmi{t,t')Hl - [ dt' [eiji-2A{t~t')6nei,2]Djm{t,t')H'i^ (44) 
Jo Jo 

where the x dependence of £ comes about only through the mean field, H{x, t), and its spatial 
gradients, because the G-invariant transport coefficients are independent of x . 
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e^Jm [Sjmlit, f) - 2At' 5n S,m2{t, t')] [H'^ + 2At'6l2H',] 
(^ijmSjmlitj t 



[S.uit, t') - 2At' 6n Sjuit, t')] [H'l + 2At'5i2H',] 
Gij2Sjii{t, t')H'i 



V. MEAN-FIELD INDUCTION EQUATION 



Applying the shearing transformation given in equations ^ and ^ to the mean-field equa- 
tion dl]), we see that the mean-field, H{x, t), obeys 



dt 
where 



rll-f 

• + 2A6aHi = (Vx£)i + r,V^Hi (45) 



It may be verified that equation (l45l) preserves the condition V- i? = 



V-if = = Hpp + 2AtHi2 = (47) 



dXp 



We now use equations (|44l) and (|46l) to evaluate V X 



/ dt' C,mi{t,t') [Hlp + 2At6p,H'i2] 

— / (it' Djmit, t') [eipqCqjl — 2A(t — t')6lieipqeqj2] [H'i^p + 2y4i(: 6plH'i^2\ 

Jo 

Expanding eipqeqjm = (Sij 5mp — 5im Sjp), the contribution from the C term is 

{Vx£)f = f dt' [C^pi - Cpu] [H'l^ + 2AtSp,H'i,] (48) 
Jo 

Evaluating the D term is a bit more involved. Again, we begin by expanding eipqGqji = 
(Sij 5ip - 5u 5jp). Then we get 

(V X £)f = [ dt' Dpm {H'^p^ + 2At6p^H'^,.^ - 2A(t - t')6,2 [H',p^ + 2At6p^H',2^] } 
Jo 

- f dt' D,^ [H'^^^ + 2At'H',^^] (49) 

^0 
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The second integral vanishes because the factor in [ ] multiplying is zero: to see this, dif- 
ferentiate the divergence-free condition of equation (l47l) with respect to Xm- Gathering together 
equations (|48l) and (l49l) . we have 



{Vx£),= I dt'[C,^i-Cmu][H[^ + 2At5^^H[^] + 
Jo 

+ f dt' D,^ {H[^^ + 2At5,^H'^^ - 2A{t - t%2 [H[^m + '^At5,^H[^^] ] (50) 
Jo 

Thus the mean field H{x, t) satisfies the mean-field induction equation, 

rlT-f 

+ 2A5,2Hi = r]V^H, + J dt' [Q^i - C^u] [K, + 2At6miK] + 

+ /* dt' D,^ {H[^^ + 2At5,^H'^,^ - 2A{t - t')6,2 [H',^„, + 2At5,^H',^„] ] (51) 
Jo 

Equation (|5T1) gives a closed set of integro-differential equations governing the dynamics of 
the mean-field, H{x, t), valid for arbitrary values of A. Some of its important properties are: 

1. Only the part of Cimiit, t') that is antisymmetric in the indices (i, m) contributes. 

2. The Djm{t, t') terms are such that (Vx£)j involves only for i = 1 and i = 3, whereas 
(Vx£)2 depends on both H2 and Hi. This means that the mean-field induction equa- 
tions (I5TI) determining the time evolution of Hi(x, t) and H^{x, t) are closed, whereas the 
equation for H2{x,t) involves both H2{x,t) and Hi(x,t). Thus Hi(x,t) (or H3{x,t)) 
can be computed by using only the initial data Hi{x,0) (or H3(x,0)). The equation for 
H2 involves both H2 and Hi, and can then be solved. The implications for the origi- 
nal field, B(X,t), can be read off, because it is equal to H{x,t) component-wise (i.e 
Bi(X,T) = Hi(x,t)). Thus, the Djm{t,t') terms do not couple either Bi or with any 
other components, excepting themselves. In demonstrating this, we have not assumed that 
either the shear is small, or that H(x,t) is such a slow function of time that it can be pulled 
out the time integral in equation 



3. When the turbulence is non helical, Cimi{t,t') = 0, but Djm{t,t') 7^ 0. In this case, there 
is no shear-current type effect, in quasilinear theory in the limit of zero resistivity. This 
result should be compared with earlier work discussed in [14|, |6|, |7D, where there is explicit 
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coupling of B2 and Bi in the evolution equation for Bi. A generalization of equation (|5T1) 



to the case of non zero resistivity has been worked out in [|lll] . It is interesting to note that 
the corresponding generalization of Cimi that appears in this case need not vanish for non 
helical turbulence. However, it is expected to vanish in the formal limit of zero resistivity, 
consistent with our result given above. 

VI. THE INDUCTION EQUATION FOR A SLOWLY VARYING MEAN-FIELD 
A. Mean EMF 

The mean EMF given in equation (l44l) is a functional of Hi and Him- When the mean-field 
is slowly varying compared to velocity correlation times, we expect to be able to approximate 
£ as a function of Hi and Him- In this case, the mean-field induction equation would reduce 
to a set of coupled partial differential equations, instead of the more formidable set of coupled 
integro-differential equations given by (l45l) and (l50l) . Sheared coordinates are useful - perhaps 
indispensable - for calculations, but physical interpretation is simplest in the laboratory frame. 
Hence we derive an expression for the mean EMF in terms the original variables Bi and Bim- The 
result may be stated simply: 



dB 

£i = aii{T) Bi{X,t) - r]imi{r) 



aii{r) = eijm / dr' [Cjmi{T,T') + 2A{t ~ T')6iiCjm2{r,T')] 







Viml{r) = tiji I dr' [Rjm{T,T') - 2A{t - T')6m2Rjl{r,T')] (52) 

which is derived below by two different methods. 

1. Method I: use of a perturbative solution for H(x, t') 

Consider the mean-field equation (1451) when S can be considered small. We introduce an 
ordering parameter e <^ 1 and consider £ to be 0(e) . Then a perturbative solution of equation (l45l) 
in the r] limit is 
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Hi{x,t') = Hi{x,t) + 2A{t-t')6i2Hi{x,t) + 0{e) (53) 

We can also consider perturbative solutions with non zero t], but using them in equation (l44l) for 
S would not be correct, because equation (|44l) was derived in the limit r/ — t- 0. We now use 
equation (l53]) in (l44l) : 



£i{x,t) = Hiei.m [ dt' Cjmi{t,t') + 2AH,eijm [ dt' {t - t')Cjm2{t,t') 
Jo Jo 

-Himeiji [ dt'Djm{t,t') + 0{e^) (54) 
Jo 

Transform to the original field variables, using Hi = Bi and Him = Bim — 2At5miBi2, which is 
given in equation (|30l) . The C terms remain unaltered and can be seen to combine to equal auBi. 
Work out the D term using the expression for Djm given in equation (1431) : 

Him / dt' Djm = [Bim — 2At6mlBi2] / dt' [Rjm + 2At' 6m2Rjl] 

Jo Jo 

= Bim / dt' Rjm — 2ABi2 I dt'{t — t')Rji 
Jo Jo 

Using the above result, and ignoring O(e^) terms in equation (|54|) . we obtain the result stated in 

equation (|52l) . 



2. Method II: Taylor expansion of B{X' , t' = t') 

This is the standard approach, although not as short as the one given above. We express 
H{x,t') = B(X',t' = t') and Taylor expand B inside the integral in equation (|44|) . As in 
equation (|23] ). 

X = {xi ,X2 — 2Atxi , X3) ; X' = (xi , X2 — 2At'xi , 0:3) 
Writing X' = X + 2A{t — t')xie2, we Taylor-expand: 
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H'l = Hi{x,t') = Bi{X\t') = Bi{X + 2A{t-t')x^e2,t') 

= Bi{X,t) + 2A{t-t')xiBi2 - {t-t')^ + ■■■ 

We now use the mean-field induction equation ([3]) to evaluate (dBi/dt). As earlier we drop the 
contributions from (V x£) and the 77 term and get 

= 2AxiBi2 - 2A612B1 + ... (55) 

ot 

Then 

H'l = Bi{X,t) + 2A{t - t')xiBi2 - {t-t')[2AxiBi2-2A6i2B,] + ... 

= Bi + 2A{t - t')6i2Bi + ... (56) 

Note that the inhomogeneous terms proportional to xi mutually cancel. It is clear, on physi- 
cal grounds that they must, because the mean EMF given by equation (l44l) is GI, and any valid 
approximation of a GI expression must preserve this symmetry. In particular, this implies that 
transport coefficients cannot depend on Xi. We now use equation (|56l) inside the time integrals of 
(|44|) . Bi = Bi{X, t) is a function of (x, t) and can be pulled out of the integrals over t'. Work out 
the C and D terms separately: 



Sf^ — €ijm I dt Cjml{t,t)Hi 

Jo 

[ dt' Cjmi [Bi + 2A{t - t')5i2Br] 
Jo 

I dt' [CjmlBl + 2A{t — t')Cjm2Bl 

Jo 



an Bi 



To calculate the D terms, we note that H'l^ = (dHi/dxm)- Since the integral over t' is performed 
at constant x, the (d/dxm) can be pulled out of the integral: 

= --^ / dt' [e,,i-2A{t-t')6ne,,2]D,Ut,t')H'i 
(y^m Jo 
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Work out 



[eiji - 2A(t - t')5aeij2] H[ = \e,ji - 2A{t - t')Sneij2] [Bi + 2A(t - t')5i2B{\ 

Then 

= —^ijlT. f dt' Djjn{t,t') 

CfXm Jo 

The quantity 

OBi f d d \ 

= r.^ — "^AtSml-^TT Bi = Bim — 2AtSmlBl2 

can be regarded as a function of (X, t) (or equivalently (a?, t)), and we are free to take it inside the 
t' integral. When this is done and the expression for Dj^ given in equation (l43l) is used, we have 

8^ = —tiji I dt' [Bim — 2At6.mlBi2] [Rjm + '2At'Sm2Rjl] 

Jo 

= —^ijlBlm / dt' [Rjm — 2A{t — t')dm2Rjl] 

Jo 

= —VimlBlm (57) 

B. Calculation of V x 5 

We need to calculate Vx£ for the mean EMF of equation (|52l) . Work out the a and r] terms 
separately. 
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(Vx£)" = eipgOiiBip = Bipeipqegjm / dr' [Cjmi + 2A{t - T')6iiCjm2] 

Jo 

= Bira [ dr' [Cirnl + 2A{t - T')6nCim2\ 

Jo 

- Bi, [ dr' [C,u + 2A{r - r')5nC,i2] 
Jo 

= Birn I dr' {Ciml — Cmil + 2A{t — t')Sii [Cim2 " Cmi2]} 

Jo 

Note that only the part of Cimi that is antisymmetric in the indices (i, m) contributes. 



^ipgVqmlBlpm Bipfy^^ipq^gji I dT [-Rjm, 1A{t T )(5}Ti2-Rjrl] 

^0 

f dr' [R,^ - 2A{t - T')6m2Rjl] 

Jo 



(58) 



(59) 



where we have used Bu = 'V- B = 0. We note that equations (l58l) and (|59] ) can also be derived 
directly from the expression for Vx£, given in equation (|50l ). This is an interesting exercise as 
it allows us to formulate an alternate criteria on when the integral equation for B can be approxi- 
mated by differential equations. We examine such an approximation further below. 



C. Approximating the integral equation directly 

It is convenient to work with the Fourier transform of H(x,t): 

H{k,t) = J d^xH{x,t)exp{-ik-x) (60) 

We also define the vector K{k,t) = {h + 2Atk2, k2, h) and K"^ = |Kp = {h + 2Atk2f + 
+ k^: note that K ■ X = k ■ x. The magnetic field in the original variables, B{X, r), can 
be recovered by using the shearing transformation, equation Q, to write {x,t) in terms of the 
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laboratory frame coordinates (X, r): 

B{X,t) = H{x,t) 



W) 



- H{k, t) exp {ik- x) 



H{k,T)exp {iK{k,T)-X) 



(61) 



From equation (|5TI) . the Fourier transformed induction equation becomes 

f)J-f ~ r ~ 

+ 2ASi2Hi = -vK^H, + tj dt' [Cmi - C^u] [H[k^ + 2At6^,Hlk2 



dt' Dj^ \ H-kjkm + 2At5jiH%km 



dt' Djm \ 2A(t - t')6i2 H[kjkm + 2AtSjiH[k2k„ 



(62) 

Let us again simplify the integrals corresponding to the C term, say T'^ and D term, say T^, 
separately. Using the definition of K(k, t), the C term simplifies to 

= iK^{k,t) [ dt'[CM-Crmi]H[ (63) 

We now assume that the mean field is slowly varying compared to the correlation time t^. of the 
turbulence and Taylor expand Hi{k^ t') about t (this assumption can later be checked for its self- 
consistency). We get 

Hiik,t') = H,ik,t) - {t-t'f-^ + ... 



Hi{k, t) + 2A{t - t')6i2Hi -{t-t' 



dt 



+ 2A5i2Hi 



+ ... 



(64) 



where in the second line we have added and subtracted a term 2A{t — t')5i2Hi. Substituting this 
expansion in equation (|63]) . the C-term becomes 

= tKm{k, t)Hi [ dt' {Cm - Cmu + 2A{t - t')6n [Ci^2 - C^s]} 



- iKm{k,t) 



m 

dt 



+ 2A5i2Hi 



dt'{t — t')[Ciml — Cmil] 



(65) 



Now consider the D-terms. Again using the definition of K(k,t) and Djm = Rjm + 
2At' 5m2Rji, we can simplify this to 



= -K,Km [ dt' [H', - 2A{t - t')6,2H[] [Rjn. - 2A(t - t')6m2Rji] 
Jo 



(66) 
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Again assume that the mean field is slowly varying compared to the correlation time Tc of the 
turbulence and Taylor expand Hi{k, t') about t. To first order in (t — t'), we have 

dili 



dt 



\H[-2A{t-t')5aH[] = Hi - {t-t' 
Substituting this expansion in equation (|66l) gives 

Tf" = -K^K^H, [ dt' [R^^ - 2A{t - t')5m2R,i 



+ 2A5,2H^ 



+ ... 



m 

dt 



+ 2A5i2Hi 



dt' {t - t')[Rj^ - 2A{t - t')5m2R, 



(67) 



The expressions for Tf and Tf given in equations (|65]) and (|67l) can be simplified. In both 
equations, the second terms are proportional to the LHS of the induction equation (|62l) . As before 
we ignore microscopic diffusion and write 



dt 



Then equations (|65l) and (|67l) can be written as, 



iK^{k, t)Hi / dt' {C^mi - C^a + 2A{t - t')6n [C,„2 - ^^2]} 



- zK^k, t) [Tf + TP] f dt'it - t') [Qml - Crml] 

Jo 



Tl" = -KjKmHi / dt' [R,^ - 2A{t - t')5rraRji] 
Jo 

+ K,Km [T^^ + T^] [ dt' {t - t') [Rjrr. ' 2A{t - t')6m2Rjl] 

Jo 



(68) 



When these equations are added together, they result in a set of three coupled linear equations for 
the unknown quantities [Tf + Tf] , [T^ + T^] and [T^ + T^] . It is straightforward to solve 
this system of equations, but the solutions assume a form which is needlessly complicated for our 
purposes. We are interested in the limit of short velocity correlations times, Tc . In this case both 
Tf and Tf are well approximated by their respective first terms: 

= iK^{k,t)Hi [ dt'{C,^i-C^u + 2A{t-t')6a[a^2-Crm2]} 
Jo 



= -KjK^H, [ dt' [R^m - ^Ait - t')6^2Rn] 
Jo 



(69) 
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These are exactly the Fourier transforms of equation (|58] ) for (Vx£)°, and equation (|59l ) for 

(vx£)^ 

We now state the conditions under which the approximations given in equations (l69l) are valid. 
Let us define the quantitites «o and r^turb as typical values of the time integrals of the velocity 
correlators, Sjmi and Rjm, respectively (for homogeneous and isotropic turbulence, ao is of order 
the magnitude of the usual a-effect, and T^turb would be comparable to the magnitude of the usual 
turbulent diffusion coefficient). For any wavenumber, K, we can define time scales, ta = {Kao)~^ 
and trj = (K'^TltuTh)^^, associated with and r^turb- When Tc is small enough such that 

Tf. <^ ta ] AT(. ^ 1 

At^ <^ ta,t^; At,>1 (70) 

then both Tf and Tf are well approximated by their respective first terms, as given in equa- 
tions (|69l ). The time scales, ta = (Kao)^^ and = (f^^r^turb) \ depend on the spatial scale, 
K^^, which is a time-dependent quantity for k2 7^ 0; at late times, K ~ |2At/c2| and this makes 
the quantities ta and decreasing functions of time. With this fact taken into account, the in- 
equalities given in equation (fTOl) translate into upper limits on the time over which the expressions 
in equation (|69| ) serve as good approximations to Tf and Tf. 

D. Mean-field induction equation 

We gather together here the results obtained in this section. When the mean-field is slowly 
varying, it satisfies the following partial differential equation: 

where 

ttimi(T) = / dr' {Cimj - Cmij + '2A{t - T')6ji [Cim2 - Cmi2\} 

Jo 

fjjmi'r) = I [ dr' {Rjm + Rmj - 2A{T-T')[6m2Rjl + Sj2Rml]} (72) 

^ Jo 
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In the above integrals Cimj = Cimjij, t'), Rj^ = Rjmij, t') etc. Some comments: 

1 . Note that otimj is antisymmetric in the indices (i, m), whereas fjjm is symmetric in the indices 
(j,m). 

2. fijm terms do not lead to coupling of any component of B with any other component. 

VII. MEAN-FIELD DYNAMICS FOR NON HELICAL VELOCITY STATISTICS 

When the velocity fluctuations are non helical, Simj{T,T') = 0, so that both CimjiT,T') and 
amj(T) vanish (in specific models of the velocity dynamics we find that the generated velocity 
fluctuations are indeed non-helical, if the forcing is non helical even in the presence of shear). Then 
the evolution of the mean-field, (over times when the inequalities of equations |70] are satisfied), is 
determined by 

Note that f]jm depends on the nature of the stirring and will, in general, be a function of time; this 
will be the case, say, for decaying turbulence. However, for statistically stationary stirring, fijm 
will become time-independent, after an initial transient evolution. 

Equation (l73l) is inhomogeneous in the spatial coordinates so, as before, we find it convenient 
to work with the new variable, H{x, t), and transform equation (1731) to the shearing coordinates 



a TT pp. TJ 



where (see eqn. 



^ = 7^ + 2At5pi A ; = + 2At5,i (75) 

aXr, OXr, OX-) \ OXn OX 



Equation (1741) is homogeneous in x but not in t, so we take a spatial Fourier transform defined 
earlier in equation (|60l) . Then H{k, t) satisfies 

+ 2A5,2H, = - [fj,m{t)K,K^ + r]K^] (76) 
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where the vector K{k, t) = [ki + 2At ^2, ^2, h) and K"^ = |Kp = {ki + 2Atk2Y + fcf + , 
as before. It may be verified that this equation preserves the Fourier version of the divergence 
condition of equation (|47l) . namely K- H{k, t) = . The solution is 



Hiik,t) = Hiik,0)g{k,t) 



H2{k,t) = H2{k,0) - 2AtHi{k,0) g{k,t) 



H,{k,t) = Hs{k,0)g{k,t) 



(77) 



where H{k, 0) are given initial conditions satisfying k- H{k, 0) = 0, ensuring that K- H(k, t) 
. The Green's function, ^(fc, t), is zero for t < and is defined for t > by 



Q{k,t) = exp 



- ds [rijm{s)KjKrr, + TjK'^) 

Jo 



(78) 



In the integrand, Kj = kj + 2As5jik2 should be regarded as a function of k and s, and the s- 
integral performed at fixed k. Then Q{k, t) can be written as the product of a microscopic Green's 
function, t), and a turbulent Green's function, Qt{k, t) : 

g{k,t) = g^ik^t) ■ g,{k,t) 



grj{k,t) = exp 



-r][k'^t + 2Akik2 + -A^A;^ 



gt{k,t) = exp[-Qjm{t)kjkr, 



(79) 



where the time-dependent symmetric matrix Qjm{t) is given by 



Qjm{t) = / ds{r]jm{s)+2As[Sj2Vlm{s)+Sm2Vjl{s)]+4:A'^Sj2Sm2s'^r]u{s)] (80) 

Jo 

in terms of time integrals of f]jm{T), which are assumed to be known functions depending on the 
velocity correlators, RjmiT, r'), as given in equation (1721) . 

The solution in the original variables, B{X,t), can be recovered by using the shearing trans- 
formation, equation Q, to write (a;, t) in terms of the laboratory frame coordinates (X, r) (see 
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equation [6TI): 



/d'^k ~ 
— /f(fc,t)exp {ik-x) 



if(fc,r)exp(zK(fc,r)-X) (81) 



(27r)3 

Equivalently, the solution is given in component form as 

5i(X,r) = j 0-^Brik,Q)g{k,T) exp {zK{k,ry X) 

/d^k r ~ 1 
^ [B2{k, 0) - 2ArB,{k, 0)J g{k, r) exp {zK{k, r)- X) 

/d^k ~ 
— B,{k,Q)g{k,T) exp{zKik,ryX) (82) 

where we have written the initial condition, H{k, 0) = B{k, 0), with k- B(k, 0) = 0. 
Some comments: 

1. The above solution for B(X, r) is a linear superposition of shearing waves, of the form 
exp {iK(k, t)- X) = exp [i{ki + 2ATk2)Xi + iA;2-^2 + ^^s-^s]^ indexed by the triplet of 
numbers (ki,k2,ks). 

2. Whether the waves grow or decay depends on the time dependence of the Green's function, 
g{k, r) = Qriik, t) ■ Qt{k, r). The first term, is known explicitly and describes the ulti- 
mately decay of the shearing waves (on the long resistive timescale), although these could be 
transiently amplified. The second term, Qt, depends on the properties of the time-dependent 
symmetric matrix Qjm{T~)- Shearing waves can grow if Qjm{T) has at least one negative 
eigenvalue of large enough magnitude. To translate this requirement into an explicit state- 
ment on dynamo action requires developing a dynamical theory of the velocity correlators, 
RjmiT, t')' because Qjmi'T) depends on time integrals over Rjm{T, r'). 

In specific cases it is possible that the velocity dynamics is such that %m(T) becomes indepen- 
dent of r, in the long time limit (this is generic when steady forcing competes with dissipation). 
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Taking the zero of time to be after this stationary state has been reached, we can do the s integrals 
in equation (l80l) explicitly and write 



Qjm{t)kjkm = t {T]jmkjkm) + ^At^ (?7lmA;m^2) + g^^*^^ (^11^2) (83) 

We can now make further statements on the dynamo growth using equation (1831) . Note that the 
hnear shear of the form that we have adopted is likely to lead to a non-zero 7712, but is not expected 
to couple the X3 component with other components, and thus we expect 7713 = 7723 = 0. Then 



- Qjmkjkm = -t [riiikl + ri22kl + 2r]i2kik2 + rjsskl)] - 2At^ [?7iifcifc2 + rii2kl] - -At^ rjukl 

(84) 

The term linear in t will dominate at early times while the term proportional to will dominate 
eventually. Thus at early times we need one of the eigenvalues of the matrix 



V 



1112 V22 

r^aay 

to be negative for dynamo growth. These eigenvalues are 



A- 



11 



^22 j 



± 



|?7ll - V22\ 



1 + 4— 



Vl2 



1/2 



A., 



^33 



(85) 



"(^11 -V22y\ 

Nonzero values of r/12 or negative values of the diagonal elements of the turbulent diffusion tensor 
favour growth at early times. Preliminary work on simple models of velocity dynamics that we 
are exploring suggests that 7722 can become negative but fjn and 7733 remain positive; this happens 
because the turbulence is strongly affected by the background shear and the velocity correlators 
are not isotropic. Thus a non-zero k2 seems to be required for growth initially. 

At intermediate times, when the term dominates we can always choose shearing waves with 
an appropriate sign and magnitude of kik2 such that 2At^ (7/11^1^2 + "ijuk'^) is negative, and there 
is growth of the mean field. On the other hand, all shearing waves with non-zero k2 will eventually 
decay, in the long time hmit t — t- 00, if r/n > 0, as then the term is negative definite. Thus it 
seems likely that the shear dynamo can have shearing wave solutions which grow for some time 
if they have non-zero X2 dependence, but which will eventually decay. As already emphasized 
above, one needs to develop a dynamical theory of the velocity correlators, for deriving more 
explicit results on dynamo action, due to non-helical turbulence and shear. It is, in general, not 
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an easy task to make analytical progress on a dynamical theory. However, in the limit of low fluid 
Reynolds numbers, a perturbative analysis is possible and the velocity correlators can be computed 
explicitly. Such an analysis has been undertaken by Singh and Sridhar, and preliminary results 
for non-helical forcing indicate that the turbulent diffusion coefficient 7722 can indeed become 
negative. Also our conclusions are based on the differential equation approximation, which is 
valid for a finite period and thus we need to solve the integral equation for the mean field evolution 
directly, to firm up the above results. 



VIII. CONCLUSIONS 



We have studied here large-scale dynamo action due to turbulence in the presence of a linear 
shear flow. Systematic use of the shearing coordinate transformation and the Galilean invariance 
of a linear shear flow allows us to develop a quasilinear theory of the shear dynamo which, we 
emphasize, is non perturbative in the shear parameter. The result is an integro-differential equation 
for the evolution of the mean magnetic field. We showed using this equation that for non helical 
turbulence, the time evolution of the cross-shear components of the mean field do not depend 
on any other components excepting themselves. This implies that there is essentially no shear- 
current type effect in quasilinear theory in the limit of zero resistivity. Our result is valid for any 
Galilean-invariant velocity field, independent of its dynamics. 

We then derived differential equations for the mean-field evolution, by developing a systematic 
approximation to the integro-differential equation, assuming the mean field varies on time scales 
much longer than the correlation times of the turbulence. For non-helical velocity correlators, 
these equations can be solved in terms of shearing waves. These waves can grow transiently at 
early and intermediate times. However it is likely that they will eventually decay at asymptotically 
late times. More explicit statements about the behaviour of the shearing wave solutions requires 
developing a dynamical theory of velocity correlators in shear flows. It is also important to directly 
solve the integral equation for the mean field as the differential equation approximation is valid 
only for a limited period. 

Growth of large-scale magnetic fields in the presence of shear and non-helical turbulence has 
been reported in some direct numerical simulations [|UL2l]. Whether we can understand these nu- 
merical results through our quasilinear theory depends on the existence (or otherwise) of growing 
solutions to the integral equation (ISTl) for the mean field. This in turn relies on the form of the 
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velocity correlators, which will be strongly affected by shear and highly anisotropic; hence it is 
difficult to guess their tensorial forms a priori, and it is necessary to develop a dynamical theory 
of velocity correlators. We cannot discount the possibility that effects we have ignored may also 
play a role. Perhaps the initial growth of the shearing wave in the mean field, for large enough 
shear, is sustained by an effect which breaks one of our assumptions. One possibility is that he- 
licity fluxes arising due to shear, turbulence and an inhomogeneous mean magnetic field [1^, 12] 



induce a nonlinear alpha effect when the Lorentz forces become strong. Another is the possible 



presence of an incoherent alpha- shear dynamo yj. 



131] in these simulations. A third possibility is 



that if even transient growth makes non-axisymmetric mean fields strong enough, they themselves 
might drive motions which could lead to sustained dynamo action; this seems remniscent of some 



of the subcritical dynamos discussed by [|l4|] . Clearly further studies of various aspects of the 



shear dynamo, particularly incorporating velocity dynamics can only be more fruitful. 
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